4 




/ 


30 ' 


STUDIES OF UNSTEADY VISCOUS FLOWS USING A 
TWO-EQUATION MODEL OF TURBULENCE 

NASA GRANT No. NAG 3-768 

SEMI-ANNUAL STATUS REPORT 

for the period 

January 1- June 30, 1987 

Submitted to 

NASA LEWIS RESEARCH CENTER 
CLEVELAND, OHIO 


73 90 

lit 


Attn: Dennis L. Huff 

Technical Monitor 


L. N. Sankar 
Associate Professor 


Oiunn-chi Wu 

Graduate Research Assistant 


School of Aerospace Engineering 
Georgia Institute of Technology 
Atlanta, GA 30332 


(AASA-CB-181293) STUDIES CE CkSlEADY 

F1CBS DSIKG A I BO- 1 fi CIIIC- B «ODEL CF 
30BBULEMCE Seiiiacnual States f<= F crt, 1 Jan 
- 30 Jun- 1987 (Georgia I DS t. cf Tech.) 11 
f Avail: NTIS BC A02/BF AC 1 CSCL 20D G3/3 


N87-27949 

Uoclas 

0C93404 





INTRODUCTION 


The objective of this research effort is to develop a two equation model of 
turbulence, based on the turbulent kinetic energy and energy dissipation, suit- 
able for prediction of unsteady viscous flows. A second objective of this re- 
search is to compare the performance of the two equation model with simpler 
algebraic models such as the Baldwin-Lomax two layer eddy viscosity model, and a 
model by Johnson and King which accounts for upstream history of the turbulent 
kinetic energy. 


SUMMARY OF PROGRESS 

The NASA Ames version of the Johnson and King model subroutine was acquired, 
and was modified for use in the Georgia Tech 2-D Navier-Stokes code, and the NASA 
Lewis version of the ARC2D code. The usefulness of the Johnson and King model was 
evaluated by performing a number of calculations with the ARC2D-Johnson-King 
solver. In Figure 1, the surface pressure distribution and displacement thickness 
variation over a RAE 2822 airfoil at 0.75 Mach number, 3.19 degree angle of 
attack and 6.2 Million Reynolds number is shown. For comparison, a set of ex- 
perimental data and results from the ARC2D-Baldwin-Lomax solver is also shown. It 
is seen that the Johnson and King model does a slightly better job of predicting 
the shock location, and the rapid rise in the displacement thickness downstream 
of the shock. In Figure 2, the velocity profiles on the upper surface at two 
chordwise locations are shown. Again, the Johnson and King model performs 
somewhat better than the Baldwin-Lomax model. 

Calculations have also been carried out for an iced airfoil configuration at 
0.12 Mach number, 1.141 Million Reynolds number for a range of angles of attack. 
In Figure 3, the surface pressure distribution at zero angle of attack as pre- 
dicted by the two turbulence models (Baldwin-Lomax and Johnson-King) are shown. 
In Figure 4, the turbulent eddy viscosity contours and streamline contours for 
the above condition are also shown. It appears that the Baldwin-Lomax model tends 
to overpredict the turbulent eddy viscosity everywhere, which results in a 
slightly shorter separation bubble. Turbulent eddy viscosity contours for the 
same airfoil at 4 degree angle of attack, shown in Figure 5, show a similar 
behavior. 


The Johnson-King model requires more computer time to obtain a final solu- 
tion than the Baldwin-Lomax model. Further, some experience with the Johnson-King 
model is needed to know at what location over the airfoil the Johnson-King ODE 
model should be turned on. Work is in progress to gain additional experience with 
the Johnson-King model. 

The two equation turbulence model (k-e) has been coded, and incorporated 
into the Georgia Tech 2-D Navier-Stokes solver. Results for the iced airfoil 
configuration using the k-c model will be presented in the next progress report. 

Appendix I gives a brief description of the equations used in the three 
turbulence models under study. 



APPEND/ A E 

TURBULENCE MODEL ( I ) 


BALDWIN-LOMAX ALGEBRAIC EDDY VISCOSITY MODEL 
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• Turbulence velocity scale and turbulence length scale are 
determined by the algebraic relations. 

t Equilibrium assumption: turbulent shear stress depends only 
on local properties of the mean flow. 

• Work well for attached and slightly separated flows. 

Some ambiguity to determine Fmax for strong separated flows 

• Easy to apply; least computational effort. 



TURBULENCE MODEL (II) 


JOHNSON-KING ODE MODEL 
- An assumed eddy viscosity distribution 
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- An ODE for the maximum Reynolds shear stress 
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Employ an ODE relating to the turbulence velocity scale. 

Non -equilibrium model: turbulent shear stress is influenced 
by " history" effects of streamwise pressure gradient. 

Perform* better for separated flow with severe adverse 
pressure gradient than the algebraic models. 

Not so robust as the algebraic models; needs more 
computational effort than the algebraic models. 



GORSKI'S k-e TWO EQUATIONS MODEL 


k = 2 <u' 2+v ' 2 > e = J/[Ou'/3y)2 + (^v'/ax) 2 ] 

kinetic energy dissipation rate 

- Outside the viscous sublayer 
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convection difussion source 

term term term 



P = production = p^C ( i y u) 2 + (^ x v ) 2 + 2 3 y u^v ) ] 


p. = eddy viscosity = C k 2 /e 
*c p 

- Within the viscous sublayer, k, e and p^ are determined 
through the algebraic relations. u 


• Employ? two PDE's for relating the turbulence length scale 
and turbulence velocity scale. 

• Non -equilibrium model; employ&more physics than the other 
two models. 

• Generally it predicts better mean flow properties and 
turbulent shear stress for separated flows than the 
algebraic model and the one-equation model. 

• Not easy to apply; requires most computational effort among 
the eddy viscosity models. 


DSTAR/C "CP 



I c I I I I I I 

1 *^.0 0.2 0 .% 0.6 0.8 1.0 

X/C 

PRESSURE DISTRIBUTIONS FOR 
RAE 2822 AIRFOIL. MINF=0.750» 
A = 3 • 1 9 » RE = 6 . 2 0 MILLION. 



L MODEL 
K MODEL 


16.1 


UPPER SURFACE OF RAE 2822 AIRFOIL 
A = 3 • 1 9 » M I NF = 0 . 7 50 * RE = 6.20 MILLION 


0/A co 0/A 


FIG. 



VELOCITY PROFILE AT X/C=0.750» 

UPPER SURFACE OF RAE 2822 AIRFOIL. 
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B-L model 



stream function contour J-K model 
j. M*0.12 alpha- O'* Re* 1.41 million 







B-L model 



stream function contour j-k model 
F<<j^ k fT M=0. 1 2 alpha»4° Re«1.41 million. 




